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I. INTRODUCTION 

The ground state of QCD at large densities is color superconducting quark matter P, 0, 0, H ■ When the effects of 
quark masses can be neglected, three flavor quark matter will be in a particularly symmetric phase called the Color- 
Flavor Locked (CFL) phase, with BCS pairing involving all nine quarks In this phase the fermion excitation 
spectrum has a gap A, and model calculations indicate that A ss 10 — 100 MeV for quark chemical potential fi « 
300-500 MeV. This phase breaks the S£/(3) co i or x SU(3) L x SU(3) R x U(1) B symmetry of QCD down to the global 
diagonal SU (3) symmetry. The lightest excitations are an octet of pseudo-Goldstone bosons and a true Goldstone 
boson associated with the breaking of the global U(1)b symmetry. At some densities, however, the strange quark 
mass may induce an appreciable stress on the symmetric CFL state, and less symmetric phases may be possible. One 
possibility is the CFLK phase, which exhibits a Bose condensate of K° in addition to the diquark condensate of the 
CFL phase-this phase breaks hypercharge and isospin symmetries 0,0]- Another possibility is the LOFF phase, which 
exhibits crystalline color superconductivity-the diquark condensate varies periodically in space, breaking translation 
and rotation symmetries |a M> 03 LLjl LUJ E3. LL-J LL3 • Most recently, a superconducting phase of three-flavor quark 
matter with non-trivial gapless fermionic excitations has been suggested fig . 

The densities at which color superconducting quark matter exists could be attained in compact "neutron" stars 
or core-collapse supernovae. It is therefore important to explore the impact of color superconductivity on observable 
aspects of these astrophysical phenomena. Investigations to that end have included studies of magnetic properties of 
neutron stars an d the equation of state of dense matter [LsL ll^. I20L l2ll 12^ . l23| . The interactions of neutrinos with 
superconducting quark matter has also been explored |24l l25l l26| . The emission of neutrinos from CFL during the 
long-term cooling epoch was studied in Ref. [27^. During this epoch the temperature is T < 10 10 K, and most of the 
excitations of CFL have small number densities. As a result, neutrino emission is highly suppressed. Complementing 
Ref. |27| is Ref. |28j . which studied the emission of neutrinos from CFL in a young, proto- neutron star. During this 
epoch, the temperature is T w 10 11 K, and the number densities of the excitations are significant. The authors of 
Ref. [2£| only studied the interactions of neutrinos with the Goldstone bosons, using the effective field theory relevant 
at energies small compared to the gap. But one might expect the fermionic excitations to become relevant at T w 10 11 
K, since these temperatures are close to the critical temperature of CFL, and the gap will be suppressed. 

These are the circumstances we analyze in this paper. We study neutral-current neutrino scattering in quark 
matter, above and below the critical temperature for CFL. We do this by calculating the quark polarization tensor, 
including the effects of pair-breaking, fermionic excitations. We also include the effects of the collective, bosonic 
excitation associated with the breaking of U(1)b by using the Random Phase Approximation (RPA) to build this 
mode out of microscopic quark-quark interactions. This excitation has the largest contribution of any of the bosonic 
modes. We begin setting up the calculation in Section II. But before proceeding to study CFL, we consider the case 
of non-relativistic fermionic modes, to make contact with previous work and to motivate the use of RPA, arguing that 
consistency with current conversation requires the inclusion of the collective mode wc do this in Section III, where we 
also go on the calculate the differential cross section for neutrino scattering in superfluid neutron matter. In Section 
IV we compute the medium polarization tensor for a one component relativistic superfluid. In Section V we compute 
the medium polarization tensor for CFL, and go on to calculate the differential cross section. We conclude in Section 
VI with reflections on this work. Further details of the calculation for the one-component, relativistic superfluid can 
be found in Appendix A. Further details of the calculation for CFL can be found in Appendix B. 
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II. PRELIMINARIES 



The goal of this article is the calculation of the differential cross section for neutral-current scattering of neutrinos in 
a dense quark medium inside a proto- neutron star. Neutrinos in a proto- neutron star have typical energies E v <C M/, 
so we may write the neutrino coupling to quarks as the four-Fermion effective lagrangian 

C z = ^^j„, (1) 

where — v 7^(1 — 75) v is the neutrino current, and — q 7 M (cy — 0^75) q is the quark weak neutral current. 
The differential cross section per unit volume in quark matter can be expressed in terms of the quark current-current 
correlation function, also called the polarization tensor, IL,,, |29ll30| : 



f d 3 a _ G% E' v [l-n v (E' v )\ 



Vd 2 n 3 dEl 32tt 3 E v [l-exp(-g /T)] 



Im (L^IV) , (2) 



where E v (E'^) is the incoming (outgoing) neutrino energy, qo = E u — E' v is the energy transferred to the medium, and 
T is the temperature of the medium. The factor [1 — exp(— qo/T)]" 1 ensures detailed balance. The factor [1 — n v {E' v )), 
where 

= e( ^)/T + r (3) 

enforces the blocking of final states for the outgoing neutrino. The neutrino tensor L^ v is given by 

= 8[2k»k v + (k ■ q)g^ - [k^q v + q^V) T ie^ afi k a q } , (4) 

where the incoming four- momentum is and the momentum transferred to the medium is q^. The minus (plus) 
sign on the final term applies to neutrino (anti-neutrino) scattering. The response of the medium is characterized by 
the polarization tensor, IP". In the case of free quarks, 

n^(Q) =-iJ |^Tr [So(P)T»So(P + Q)n , (5) 

where Sq(P) is the free quark propagator at finite temperature and chemical potential, and r M = 7 M (cy — 0,475). The 
inner trace is over Dirac, flavor, and color indices. The free quark propagator is [31| 

So( P )=% /(g±||A> + g^A pT „), (6) 

where c, d and /, g denote color and flavor indices, respectively, and 



A± = -(l±7o7-p) (7) 



are the positive and negative energy projection operators. The energies in the propagator are measured relative to 
the Fermi surface. To wit, E p = p — ji for (massless) particle states, and E p = p + /1 for anti-particle states. 

Since we are interested in the interaction of neutrinos with a superconducting medium, we introduce the Nambu- 
Gor'kov formalism |32t l33j| . which allows the incorporation of diquark correlations into the quark propagator. This is 
done by artificially doubling the quark degrees of freedom by introducing charge conjugate field operators qc and q~c, 
defined by 

q = Cq~ c T and q = q^C, 

where C = 17072 is the charge conjugation matrix (with C = —C^ 1 = —C T = — C 1 * and C^C^ 1 = — 7J). The 
Nambu-Gor'kov field is given by 

* = ( ^ ) and *=(q qc). (8) 
In terms of these fields, the weak interaction Lagrangian, Eq. becomes 

£ z = ^ 7/1 (l-75>*r£* (9) 
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where the neutrino-quark vertex in the Nambu-Gor'kov space is 

z ~\ -^{c v + c Al5 ) J ■ [W) 

The Nambu-Gor'kov propagator that includes BCS diquark correlations in the mean-field approximation is given by 

s(p) = ( £S %X\ ) , (id 



where 



G-(P)J ' 



^0 ?p /*0 ?p 



S+ W = ^2 7 5 A p , E-(P) = 75 A+ , (1.2) 

Po s P Po ?p 



and 



We have suppressed any color-flavor indices, and we are neglecting the contribution of antiparticles. 

III. A NON-RELATIVISTIC DETOUR 

Although we will ultimately be interested in relativistic, superconducting quark matter, we digress to consider a 
one component system of non-relativistic superfluid baryons, such as neutron matter with 1 Sq pairing, which is itself 
relevant to neutrino transport in neutron stars and supernovae. This will allow us to make contact with the vast 
body of published work on superconductivity in non-rclativistic systems and will serve as a pedagogical prelude to 
the relativistic case. 

For non-relativistic fermions the structure of the Nambu-Gor'kov propagator greatly simplifies and is given by a 
2x2 matrix 36]. We shall consider a simple model for the pairing interaction. In this case the non-relativistic 
Hamiltonian, in terms of the fermion creation and annihilation operators, is 

H - [J.N = ^ E p a l a p + G^ata^pa-p-ap, (13) 
p p 

where N = J] a^a is the number operator, E p = p 2 /2M — /i, M is the mass and /i is the non-relativistic chemical 
potential. In the mean field or BCS approximation, we replace the pair operator by the space and time independent 
classical expectation value. This defines the constant field A* = — GJ2 P a ^-p* an d A = — GJ2 P apa-p- The solution 
of the gap equation, namely d(H — /iN)/dA = determines the magnitude of A. The propagator is given by 

W = ( P(p] gV) ) < (14) 



where 



G+(P) = ^±|i, G-(P) = *-=* 
Pa ~ £p Po ~ £p 



Po - £p Po - £p 



The quasi-particle energy is given by 



P 2 



£ P = \/El + A^ , where £ p = ^ -/x. 




FIG. 1: Feynman diagram corresponding to the vertex equation (EqEU- The squares represent the dressed weak interaction 
vertex, and dashed lines represent the interaction between the fermions. The solid, directed lines represent Nambu-Gorkov 
propagators. 

The neutrinos couple to the weak neutral-current of the neutrons, which has a contribution from both the vector 
current and the axial vector baryon current (see Eq. Q). The neutral-current carried by the neutrons is given by 

in = % * f (p) (c v j ll (p + q,p) + c A rf(p + q,p)) %(p + q), 

where cy = 1/2 and ca = 1.23/2. In the non-relativistic limit 7/i (p + q,p) — {t 3 , 0) and 7/f (p + q,p) = (0, a 1), with 
t 3 = diag(l, —1) and 1 = diag(l, 1) in Nambu-Gor'kov space, and a is the Pauli matrix in spin space. To begin, we 
focus on the vector current, though we will return to the axial vector current. 

Conservation of the vector baryon current means that the vector response function 

rv /„ _ . f d4 P 



n;,(go, q) = ~% f Tr [S(p) 7|i (p,p + q) S(p + q) 7 „(p + q,p)] (16) 



must satisfy the conservation equation 



rV = q v = . (17) 



The polarization function in Eq. I|16|) does not satisfy this equation. This violation can be traced to the use of the 
dressed mean-field Nambu-Gor'kov propagator on the one hand, and the use of the bare vertex on the other. To ensure 
local baryon conservation, we must use a dressed vertex, T(p-\-q,p), which must satisfy the generalized Ward identity 
for the superfiuid: 

q^f^p + q,p) = TzS-Hp) - S-\p + q)r 3 . (18) 

That the bare vertex 7 M (p + q,p) is not sufficient to satisfy current conservation was realized shortly after the devel- 
opment of the theory of superconductivity by Bardeen, Cooper and Schrieffer |3?| in independent articles by Bardeen 
[38j , Bogoliubov [33, |4(j , Anderson [41j , and Nambu [33| • It is this correction that naturally leads us to incorporate 
the Goldstone mode. 

The dressed vertex T compatible with Eq. I|18|l satisfies the following integral equation 

Tn(p + q,p)=^(p + q,p) + Gi J 0^ r 3 S(p' + q)t^' + q,p')S{p')T 3 , (19) 

where G is the 2-body contact interaction defined by the Hamiltonian in Eq. (jl3JI . The behavior of the dressed vertex 
function in the long wavelength limit can be inferred from the Ward identity in Eq. jTSj. The right-hand side of 
Eq. I|18|) . in the q — > limit, reduces to 

lim (t 3 S-\p) -S-\p + q)r 3 ) = 2i Ar 2 . (20) 

The Ward identity therefore dictates that the vertex function be singular for zero energy transfer (qo = 0) when 
the three momentum transfer vanishes (q — * 0). This singularity is characteristic of the presence of a zero energy 
collective mode - the Goldstone mode expected on general grounds O ne ma y obtain the dispersion relation 

for the Goldstone mode by inspecting the pole structure of the solutions of Eq. (|19|) . For the simplified interaction 
considered here, the integral equation for the vertex function is represented by the Feynman diagrams shown in Fig.^ 



Calculating the vertex function by retaining only the terms shown in the figure is equivalent to the Random Phase 
Approximation (RPA) |3(j. Note that it contains both the direct and exchange diagrams (the second and third terms, 
respectively, on the right-hand side of Fig.^J. We only need to solve for Tq, rather than all four components of T^, 
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since for small q the dominant response of non-relativistic neutrons is due to the coupling of the baryon density to 
neutrinos; the coupling of the velocity-dependent components of the baryon current to neutrinos is suppressed by a 
factor q/M and may be neglected. An analytic solution for Tq exists and is given by 

To = — t 3 + - 7H with j H = I 

Xo X \ iV 2 u 

where 





= 2(1111:12 + 


1111:121122:11 — 


IIl2:XXnx2:X2) > 


N 2 


= 2(IIl 2: ll + 


1112:111111:22 — 


niX:X2lIx2:X2) > 


Xo 


= (1 +11x2:12 


-n U:U ), 




V 


= i - n 12: i2 - 


"1" rill:22 + II22: 


11 + 1111:221122:11 


X 


= 4n U: x2nx2 


:liril2:12 ~ (1 " 


" IIi2:X2) n i2: 12 



- (1 + ^1:22) (-1 + 2n?2:XX + ni2:l 2 ) + n 22 :xx 

- (riii:22 ( — 1 + IIx2:X2) + IIx2:X2) Il22:Xl 

- 2n? 1:X2 (i + n 22: ii) + n U:U (1 - n 2 12:12 + u 22:11 + n U:22 (i + n 22: ii)) 

H-ijM = GiJ -^-^ Sijip + q) Skiip), 

and S(p) is the mean field propagator appearing in Eq. 1)14(1, The complex function x{Q0iQ) — when qo = q c s , 
where c s = pp / (VSM) is the sound velocity of the neutron gas (for q Pf)- The vertex is singular at this point, 
indicating the presence of the Goldstone mode in the response. Note that the singularity does not appear in the vertex 
proportional to T3 in NG space. The off-diagonal nature of the singularity is characteristic of fluctuations of the phase 
of A - which is the Goldstone excitation. Further, at small energy transfer qo <ti A the expression for \ simplifies and 
we find that x — Xo'D- The Goldstone singularity occurs when T> = and xo 7^ 1 is a result of additional screening 
corrections included in RPA in both the quasi-free and Goldstone-mode responses. 
In RPA, the polarization tensor which satisfies the Ward identity is 

n£ PA (<?o, q) = J Tr \S{p) f S(p + q) r 3 ] (21) 

The simple form of the solution for the dressed vertex for go <SC £f allows us to write 

n^- 4 (q , q) = — < (q , q) + H H (q , q) , (22) 
Xo 



where 



and 



U'{ r ^,u)= -' ! Tr [S (p) t 3 S(p + q)n] (23) 

n H (qo,q)= J -0^Tr[S(p)r 3 S(p + q) lH }. (24) 

U q f is just the quasi-free or mean-field response, and n# is the response due to the Goldstone mode. We plot -Irally 
and -Imll^" 4 in Fig. [3 as the dashed line and the solid line, respectively, at various temperatures and ambient 
conditions appropriate for the neutron superfluid in neutron stars, namely, for chemical potential ^ — p 2 F /2M = 30 
MeV, and gap A = 5 MeV, where M and pp are the mass and Fermi momentum, respectively, of neutrons. Note that 
the quasi-free response has support only for qo ^ 2A at T < A, due to the gap in the fermion excitation spectrum. 
Eq. (1221 has xo 7^ 1 because correlations result in screening of the coupling of quasi-free excitations to the external 
current. In the weak-coupling limit (A fi), this correction is small, but it becomes important even at moderate 
coupling strength, as seen in the figure. At small but non-zero temperature (T <C A), Im T>(q ,q) ^ ~ exp(— T/A) 
due to a small contribution from thermal quasi-particle excitations. This reduces the singular behavior to a resonant 
response (Lorentzian) with a small but finite width. With increasing temperature, the collective mode gets damped 
as the single-pair excitations become increasingly important. This trend is clearly seen in the Fig- EI Finally at T c , 
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FIG. 2: We plot -Imll as a function of the dimensionless variable qo/(qc s ), where c s = pjr/(v3M) and with fixed momentum 
transfer q — O.lpF- The dashed line is the quasi-free vector response, -Imll^/, and the solid line is the full RPA (that is, quasi- 
free plus collective with screening corrections) vector response, -ImII^ PA . The spike at qo/{qc 3 ) ~ 1 is the Goldstone mode 
response. The dot-dashed line is the RPA axial response, -ImIl2 PA . We use neutron chemical potential /i = p 2 F /(2M) — 30 
MeV. The gap at zero temperature is A = 5 MeV, and the critical temperature is T c — 2.8 MeV. 



both the pair-breaking and Goldstone mode excitations disappear and the single-pair response dominates. The slight 
enhancement seen at small go is a characteristic of attractive correlations induced by the strong interactions in the 
normal phase |30|. 

The response in the axial vector channel can also be calculated using the method described above. In this case, the 
vertex does not exhibit any singular behavior, since there is no Goldstone mode in this channel. In a homogeneous 
and isotropic system, the axial polarization tensor is diagonal, with equal components, 11^ = n^ 2 = II^ 3 = 11^. An 
explicit calculation shows that the axial polarization tensor in RPA is 

U 1 l PA (qo,q)= -H"J(q ,q) (25) 
Xo 

where 

Tl g A f (Qo,q)= -» J ^Tr[S(p)iS(p + q)i}. (26) 

The dot-dashed lines in Fig. [21 show the behavior of the imaginary part of the II PPj4 (go , q) at fixed q as a function 
of go- The differential cross section for neutrino scattering in the neutron matter can be written in terms of the 
imaginary part of the polarization tensor. The differential cross section per unit volume for a neutrino with energy 
E v to a state with energy E' u — E v — qo and scattering angle is given by 

d 2 a Gf , 2 l-n v (E' v ) 

V d cos 9 dE' v 4tt 3 " 1 - exp(-g /T) 

x [4 (1+COS0) lmIi^ PA {q ,q)+c 2 A (3-cosfl) Im II« PA (g , g)] . (27) 
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FIG. 3: We plot the differential cross section for neutral-current neutrino scattering in superfluid neutron matter as a function 
of the dimensionless variable qo/(qc s ), where c s = pf /(y3M) ~ 0.15. The solid curves are the RPA results which include the 
Goldstone contribution and the dashed lines are the mean-field results. The incoming neutrino energy and the momentum 
transfer q were set equal to -n x T - typical of thermal neutrinos. The spike at qo/(qc a ) ~ 1 is the Goldstone mode response. 
We use neutron chemical potential /i = p%/(2M) = 30 MeV. The gap at zero temperature is A = 5 MeV, and the critical 
temperature is T c = 2.8 MeV. 

These expressions permit calculation of the neutrino opacity in superfluid neutron matter at arbitrary temperature. 
The results are shown in Fig. El The incoming neutrino energy E v and the momentum transfer q were set equal to ttT. 
This is typical for thermal neutrinos. Note that the Goldstone mode continues to play a role even when T ~ T c /2. 
Note also that with increasing temperature the gap equation yields smaller gaps - seen in the decreasing threshold for 
the quasi-free response. Although it is easy to deduce the relative importance of RPA corrections to the response from 
Fig. we present a table to quantify the differences. The table below provides a comparison between the differential 
cross sections integrated over the energy transfer (area under the curves shown in Fig. [3J • 



T 


da qf /Vdcos0 x 10" 5 m" 1 


do-RPA/VdcosO x 10 5 m 1 


0.5 T c 


0.1 


0.5 


0.85 T c 


2.6 


4.9 


1.1 T c 


12.8 


12.9 



Although scattering kinematics do not probe the response in the time-like region, neutrino pair production does 
arise from time-like fluctuations. As a result these same polarization functions can be used in calculations of the 
neutrino emissivities in superfluid matter - a process that is commonly referred to as the pair-breaking process |13| . 

At low temperature and low energy, the response is dominated by the Goldstone mode. In this regime it is 
appropriate to use the low energy effective theory involving only the Goldstone mode. The effective Lagrangian for 
the U(l) mode is given by 

Cgb = f -f(d Q Ud Q U^ - c 2 s diUditf) (28) 
o 

where U — exp(2iH/ fa), H is the Goldstone field, and fjj = M Pf/tt 2 is a low energy constant which is the equivalent 
of the pion decay constant in the chiral lagrangian. We can compute the coupling of the H mode to the neutrinos by 
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matching to the weak current in the microscopic theory. We find that the amplitude for the process v — > H v is given 
by 

A v ^ Hv = ^= c v }h d H vfv . (29) 

Using Eg. l[29l). we can compute the differential cross section for neutrino scattering due to the "Cerenkov" process 
v^Hv |2j|. We find that 

G| /A 4 K 2 , X ~ n r (K L (l + cos^) go %o - c s g) • (30) 



V d cos 6 dE' v An 2 JH v v 1 - exp(-g /T) 

The low energy and low temperature limit of the RPA response should agree with the above result obtained using 
the effective theory. We show that this indeed the case. We begin by noting that the RPA vertex when go "C A and 
q <C kp can be written as |3^| 

- 2A go . 

Substituting this result in Eq. IL'l'l) we find for qo > 

Im n« PA (q Q , q) = ~^q S(q Q - c s q) . (32) 

The quasi-free response is exponentially suppressed at low temperature, and using Eq. (|32(l it is easily verified that 
the differential cross section in RPA agrees with Eq. 1|3(JI) . 



IV. RESPONSE OF A ONE-COMPONENT, RELATIVISTIC SUPERFLUID 

In the last section we argued that a consistent calculation of the medium polarization tensor must include not only 
the quasi-free response of the medium associated with the pair-breaking excitations, but also the collective response 
associated with the massless excitation H arising from spontaneous breakdown of U(1)b- We then calculated the 
medium polarization tensor for a one-component, non-relativistic superfluid - superfluid neutron matter. In this 
section we will calculate the medium polarization tensor for a one-component ultra-relativistic superfluid - as a warm 
up for the CFL case, and also for its own sake - including both the quasi-free response and the collective response 
associated with the Goldstone mode. 

The interactions of interest are 

Ant = £z + £s- (33) 

Here, Lz represents the interactions of neutrinos with the medium and is given by Eqs. H9I1()(I . For concreteness we use 
Cy = — 1 and ca = —1.23. Cs represents the strong interactions in the medium that give rise to the superfluid mode 
H associated with the breaking of U{\)b- Since we expect these interactions to have the form H^qij^qc + qci"(5q) 
(see Refs. ^HES0|)i we use a four-quark interaction (qij^qc + <Zc*75'z) 2 ■ To express this in terms of Nambu-Gor'kov 
spinors, we define 

iW° *»y (34) 

Then the relevant quark self-interactions are 

C s = G (4T H *) 2 . (35) 

We will calculate two contributions to the medium polarization tensor, depicted in Fig. 01 using the Nambu-Gorkov 
propagator given in Eq. I|llfl . Our formalism, notation and conventions closely follow those of Ref. [3l| . 

The first term on the right-hand side corresponds to the quasi-free response of the medium. The second term 
captures the collective response of the system associated with the superfluid excitation H . The collective response 
can be expressed as the RPA sum depicted in Fig.[SJ We will first calculate the quasi-free response, then the collective 
response. 



FIG. 4: The two contributions to the medium polarization tensor. The first term on the right-hand side of Fig. 0] corresponds 
to the quasi-free response of the medium to the neutrino probe. The second term on the right-hand side corresponds to the 
collective response associated with the Goldstone mode H. 

FIG. 5: The collective response of the medium can be evaluated in the Random Phase Approximation by summing a series of 
diagrams. 



'/ 10? [ A ( fi p. £ p+<.) {(4 + 4)T(p. p + <i) " 2 ^ c ^f (p, p + q)} 



The quasi-free response of the medium comes from first term on the right-hand side of Fig. 0] This diagram makes 
the following contribution to the polarization tensor: 

n£(Q) = 

1 r d 3 P 

2 J (2^ 

+A(-E p , -£ p+q ) {(4 + c 2 A )Tr(p, P + q) + 2c v c A -Wl v (p, p + q)} 
-A 2 B(E p , £ p+q ) { (c 2 y + (p, p + q) + (p, P + q)] 

+ 2 CyCA [Vr(p,P + q)-V^(p,p + q)]}]. (36) 

The Matsubara sum led to the quantities 
A(E',E) = 

■(t;-E)(?-E>) (Z + E)(? + E'y 



+[1 _„ (?) _„ (f)U '(i±M_(i_M )) 



and 



B(E',E) = 
1 



7v{W0-"(O] 



(2£)(2£') LL v " VWJ V?o+e-e + 



with 



? = v 7 ^' 2 + A 2 and £ = v 7 ^ 2 + A 2 . 
The trace over Dirac indices led to the quantities 



1 + Pi ■ P~2 if /i = and ^ = 

,,, , ±(pi + K) v if /J = and , . 

*± iPijP2j - \ («jyj 

±(pi+p 2 ) p ii A* 7^ an d ^ = 

5^(1 - K • pS) + (K^CK)" + (^"(Pi)" if I" ^ and i/ ^ 



and 



Tie^ 0a (pl - pS) a if MT^O and z/ ^ 

and 



(38) 



wr(Pi, P2 )= ie Tl ( S a( ^ !! T'T" (40) 



"± (,Pl)P2j - < M „, . (41) 

-7£ pi,p 2 ) if v ^ 
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o 



FIG. 6: The diagrams of the Random Phase Approximation are a geometric series. 



and 



Vf (pi, P2 ) = 



-W^(pi,pa) for v = 
W^(pi,p 2 ) for v^O 



(42) 



Details of this calculation are given in Appendix A.l. 

Now we want to calculate the contribution to the polarization tensor from the superfluid mode H . This contribution 
is the sum of a series of terms, expressed diagrammatically in Fig. [S] The series is a geometric series, and its sum is 
expressed diagrammatically in Fig. The contribution of this geometric series to the polarization tensor is 



where 



G I»(Q)I"(-Q) 
2 g(Q) 



iA x 



('v 



d 3 p 
(2^)3 

C(E P , £ p+q ) Uf(p, p + q) + C(-E p , -£ p+q ) l£°(p, p + q) 



(43) 



(44) 



with 



C(E',E) 



E-E'-Z + g E-E' + £-£' 



(20(2f) 

+[l-„(0-n(O](- 



E - E' + £ + £' E — E' — £ — 



}■ 



(45) 



and 



ff (Q) = 1 + 2G 



d 3 p 
(2^)3 



(1 + p-p + q) 



1 



(2£ p+q )(2£ P ) 



{[n(£ p+q ) - n(Cp)] (£ P+q £ P ~ £ P+q £p - A 2 ) 



1 



1 



v go + Cp+q — go - £p+q + £p 

[1 - ™(£p+q) - "Kp)] fe> + q£p + ^P+q^P + A 2 

1 1 
go - Cp+q — £p go + Cp+q + £p 



(46) 



Details of the evaluation are given in Appendix A. 2. We should check that g(Q) vanishes for qo = and q = 0; this 
will indicate the presence of a massless excitation, namely, the H boson. Setting q$ = q = 0, we need 

= l-2 G /|^[l-2^ p )]l. (47) 

This easily verified by noting it is the gap equation obtained by minimizing the free energy: dft/dA = 0. 

In Fig.|7|we plot -Im II 00 for various temperatures as a function of qa/(qc s ), where c s = 1/V3, which is the expected 
velocity of the Goldstone mode. The dashed line is the quasi-free response, and the solid line is the full (that is, quasi- 
free plus collective) response. The momentum transferred to the medium is fixed at q = 60 MeV. We use chemical 
potential fi = 400 MeV. The gap at zero temperature is A = 50 MeV, and the critical temperature is T c — 28.3 MeV. 
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FIG. 7: We plot -ImII 00 for a one-component, relativistic superfluid at various temperatures, as a function of the dimensionless 
variable qo/(qc s ), where c B = l/VS, with fixed q = 60 MeV. The dashed line is the quasi-free response, and the solid line is 
the full (that is, quasi-free plus collective) response. We use chemical potential fi = 400 MeV. The gap at zero temperature is 
A = 50 MeV, and the critical temperature is T c = 28.3 MeV. 



The upper left panel illustrates T/T c = 0.50. At this temperature the gap is A = 47.9 MeV. We find a narrow peak 
about qo/(qc s ) = 1, confirming our expectation about the Goldstone mode. The threshold seen at qo/{qc s ) — 2.76 
corresponds to qo = 95.8 MeV = 2A. The middle panel illustrates T/T c — 0.85. At this temperature the gap is 
A = 32.2 MeV. The peak associated with the Goldstone mode has widened, and its height has shrunk relative to the 
response at qo = 2A. Also note the nonzero response for qo/(qc s ) < 1.73, corresponding to go < ?■ The rightmost 
panel illustrates T/T c = 1.10. At this temperature the gap is A = 0. Only the qo < q response is seen, as expected. 
Having studied the medium polarization tensor for a one-component, relativistic superfluid, we turn to the case of a 
multi-component, relativistic superfluid - color-flavor locked quark matter (CFL). 



V. RESPONSE OF CFL 



We finally study CFL. The interactions can still be written as Eq. (|33|l . except that for three-flavor quark matter 
we should use 



/ 1 4 • 2 n 
I - — - sm % 



cy = 



2 3 



V 



o 

- \ + | sin 2 9 W 




\ 



(48) 



\ + fsin 2 ^ J 



where the weak mixing angle is sin 9\y ~ 0.231, and 



ca 



(\ 
-| 

V° "3 



(49) 
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Also, the four-quark interaction that gives rise to the H in CFL is 

C s = G(*r ff *) 2 , (50) 

using 

r H = ( . ° ** M ) with (M)f g = eV, (51) 

where c, (i are color indices, and /, g are flavor indices. In CFL the nine quarks form an octet with gap As = A and a 
singlet with gap Ai = 2A. We are neglecting condensation in the color 6 channel, keeping only the condensate in the 
3 channel. The propagator involves the energies £ m p = \/{p— /i) 2 + A^. The propagator also involves the following 
color-flavor matrices: 

(Pi)f. = (Pi)f g = \r& (Ps)f a =S fg S^- 1 -5}5 d g , and (P 8 )f g = -Sfa + ^ 
Now we can write down the quark propagator: 

S(P) 

where 



G+(P) S"(P) 

a+(p) g-(p) 



G+ ( P ) = { ZJ^T- + -2 P % } (Po + £ P )A+7o, 




G ~( P ) = < Zo^F^ + ^ PS , 2 > (Po - P P )A p7 o, 



-+, i A|P, A 8 P 8 ^ 

~ ( P ) = - i „2_,2 + _2 .2 ? 7 5 A p ; 



and 



— i AiPi A 8 P 8 i 

a (p) = tft^ + ft^ 75V 



We have neglected the contribution from antiquarks. 

We now calculate the contributions to the medium polarization tensor from the quasi-free response and the collective 
response, the two terms depicted on the right-hand side of Fig. 0] 

The quasi-free contribution to the polarization tensor comes from the first term on the right-hand side of Fig. 0] 
and its value is 



n ^ (0) = \ fj^fl £ (A(i?p,£p +q ,A m ,A„){i?Wrf(p,p + q)-€wr(p,p + q)} 

J ^ ' m,ne{l,8} 

+A(-E p , -P p+q , A m , A„) {p« Tr (p, P + q) + BZl MT(P, P + q)} 
-A m A n P(P p , P p+q , A m , A„) {p« [Wf (p, p + q) +U^(p, p + q)] 

+ P< 2 > [Vf (p, p + q) - V^( P , p + q)]}) . 

(52) 

The Matsubara sum led to the quantities A and P, where 
A(E',£,A',A) = 

1 ^Ho-n,.a' (( - E)(( '- E ' } «+ £ >«' +£ - 



(20(2C') LLV y v /J V + <?o-£ + £' 

+[1 _„ (0 _„ (O1 ((i±M_(i_M l) 
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and 



B(E' , E, A'. A) = 

(2^r) {[n(0_n(O \*+€-€' 

+ [1 - n(0 - n(C)] f ^7 + TrW ) : 



(54) 



with 



f ' = V(^') 2 + ( A ') 2 and £ = V^ 2 + A 2 . 
The trace over color-flavor indices led to the quantities 

^(i) _ / 0.0556 0.2865 \ r{2 ) _ ( 0.0556 0.2391 \ 
mn ~ I 0.2865 2.4502 J ' m ™ ~~ y 0.2391 2.1182 J ' 

^) = ( 0-0556 -0.2865 \ ^ = / 0.0556 -0.2391 \ 

mn \ -0.2865 -0.1286 I \ -0.2391 -0.0338 J 

Note that these matrices include off-diagonal components: the weak interactions can scatter the singlet fermionic 
quasi-particle into one of the octet fermionic quasi-particles, and vice-versa. We discuss the color-flavor trace further 
in Appendix B.l. 

We now want to calculate the contribution to the quark polarization tensor from the collective response of the 
medium. This is depicted in Fig. El Its value is 



where 



I"(Q) = / E ^rnR^ [<5(£p,£W q) £ mp ,£ m p +q )Wf ( P ,P + q) 



d 3 p 

i 

m£{l 

' 1 '-p • q- p- ^>m p+q 



+ CK,-£ P+q ,( m p,( mp+q )^(p,P + q) , (57) 



with 



— {[n(0 - n(£ )] 



+ l-n(£)-n(£) t — ^ + TTTTi ( 58 ) 



("3") 

The quantity R rn arises from the color-flavor trace. Its value is 



Also, 



where 



i?i 3) =_I and 4 3) = -^- (59) 



3(Q) = 1 + 2G J ^(l + p-f+q)J2 D (Ep,E p+tl ,A m )R^, (60) 



D(E',E,A) = 

Pf)W) {'" (f) - " (r)1( «' - E£ ' - A>) (s+TT - ST^T?) 



i | -»'i'-»"''i^ ,+c£ '- A, )( s 4^?-s+?+? 11 - 
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and 

Rf ) = 4 and = 8. (62) 

We discuss the color-flavor traces further in Appendix B.2. 

The quantity g(Q) should vanish when qo = q = in order to describe a massless excitation. That is, we must have 

= 1-2 Of ^{[l-2n(£ lp )]-^ + [l-2n(£ 8p )]-M. (63) 

J ( 27r ) I Sip ?8pJ 

This is the gap equation. It gives a critical temperature T c w 0.71 A(T = 0), in agreement with Ref. |4Sj . 




FIG. 8: We plot -Imn 00 for CFL at various temperatures as a function of qo/{qc a ), where c s = l/v3, with fixed q = 60 MeV. 
The dashed line is the quasi-free response, and the solid line is the full (that is, quasi-free plus collective) response. We use 
quark chemical potential /j, — 400 MeV. The gap at zero temperature is A = 50 MeV, and the critical temperature is T c = 35.7 
MeV. 

In Fig. [S] we plot -Im LI 00 for various temperatures as a function of qo/{qc s ), where c s = 1/V3, which is the 
expected velocity of the Goldstone mode. The dashed line is the quasi-free response, and the solid line is the full (that 
is, quasi-free plus collective) response. The momentum transferred to the medium is fixed at q — 60 MeV. We use 
quark chemical potential fi — 400 MeV. The gap at zero temperature is A = 50 MeV, and the critical temperature 
is T c = 35.7 MeV. The leftmost panel illustrates T/T c = 0.50. At this temperature the gap is A = 46.7 MeV. As in 
the one-component case, we find a narrow peak about qo/{qc s ) = 1, confirming our expectation about the Goldstone 
mode. There is a threshold at qo/(qc 8 ) = 2.69, corresponding to q Q = 2A, and there is an additional small response 
for qo/{qc s ) < 1.73, that is, qo < q. The CFL response function is much more busy than in the one-component case 
(see Fig. 0) because of additional effects at go = A,3A,4A, arising from the fact that some fermionic quasi-particles 
in CFL matter have a gap A and some have a gap 2A. The upper right panel illustrates T/T c — 0.85. At this 
temperature the gap is A = 30.0 MeV. The peak associated with the Goldstone mode has widened, and its height 
has shrunk relative to the response at go = 2A, but it overlaps with the qo < q response, which has grown. The lower 
panel illustrates T/T c = 1.10. At this temperature the gap is A = 0. Only the qo < q response is seen, as expected. 

In Fig. El we plot the differential cross section for neutral-current neutrino scattering in CFL. The energy of the 
incoming neutrino is E v = 7rT, the typical energy of a thermal neutrino. The momentum transferred to the medium is 
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FIG. 9: We plot the differential cross section for neutral-current neutrino scattering in CFL as a function of the dimensionless 
variable qo/(qc s ), with c s = l/VS and with momentum transfer set to the typical incoming neutrino energy q — E v — irT. 
The spike at qo/(qc s ) ~ 1 is the Goldstone mode response. We use quark chemical potential fi = 400 MeV. The gap at zero 
temperature is A = 50 MeV, and the critical temperature is T c — 35.7 MeV. 

set to q — E v . Kinematics demand —1.73 < qo/(qc s ) < 1.73. The parameters are otherwise as in Fig- EI The leftmost 
panel is again T/T c = 0.50. There are prominent peaks near qo/(qc s ) = ±1 corresponding to the Goldstone mode. 
The expected feature at go > 2A, corresponding to excitation of a particle-hole pair, does not lie in the kinematically 
allowed region of q . There is a small peak, however, at qo/(qc s ) — —1.44, or — q n = A, corresponding to scattering 
of a singlet quark, whose gap is 2A, into an octet quark, whose gap is A. The middle panel illustrates T/T c = 0.85. 
The Goldstone peak has widened, and the threshold at qo > 2 A - corresponding to qo/(qc s ) > 1.09 - now falls in the 
kinematically allowed region. The rightmost panel illustrates T/T c = 1.10. This is the result for free quarks. 



VI. CONCLUSION 



We have calculated the differential cross section for neutral-current neutrino scattering in superfiuid neutron matter, 
plotted in Fig. [3] and in color- flavor locked quark matter, plotted in Fig. EI under conditions relevant to proto-neutron 
stars. Our results apply above and below the critical temperature. In both of these regimes our model for the 
interaction in the medium includes the dominant contribution to the cross section. Above the critical temperature, 
this comes from the fermionic excitations, which become the pair-breaking excitations below the critical temperature. 
Well below the critical temperature the dominant contribution comes from the massless bosonic mode associated with 
the breaking of U(l). Although we presented results results for scattering cross sections with space-like kinematics, 
our polarization functions extend into the time-like region where pair-breaking (and recombination) is the dominant 
source of the response. These could be employed in calculations of the neutrino emissivity. In particular, we have 
demonstrated the importance of vertex corrections in these region. Our results suggest that earlier calculations of 
the neutrino emissivity from the pair recominbation process in superfiuid neutron matter an d quark matter [25| , 
which ignore these vertex (RPA) corrections, need to be revised. 

The analysis presented here is based on mean-field and the random phase approximation (RPA). Its validity is 
restricted to weak coupling, A <C fi. For strong coupling we can expect the response to differ quantitatively. In 
particular, screening corrections which were discussed earlier could be significant. Another drawback is our use 
of simplified interactions to describe superfiuid neutron and quark matter. Our focus was to explore the role of 
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superfluidity and this motivated our choice of a simple zero-range s-wave interaction. In reality, the nucleon-nucleon 
and quark-quark interactions are more complex. These will induce additional correlations which will affect both the 
gap equation and the response. (For a recent review on the role of strong interaction correlations on neutrino opacities 
see Ref. ^^-) Given the non-perturbative nature of these corrections it is difficult to foresee how large they may be. 

Although our results are valid both above and below the critical temperature, we have implicitly assumed that 
the transition is a BCS-like second order transition. Several caveats must be borne in mind when using our results 
near T c . In the real system this transition may be first order either due to gauge field fluctuations [50j or stresses 
such as the strange quark mass and electric charge neutrality |5l|. Also, fluctuations of the magnitude of the order 
parameter dominate in a region called the Ginsburg region around T c . In strong coupling, the size of this region 
could be significant fraction of T c [52|. Our approach captures some of these fluctuations thr oug h RPA. Nonetheless, 
a Landau- Ginsburg approach - an effective theory for |A|, is more appropriate in this regime |53|. In particular, there 
are precursor fluctuations just above T c , which are not included in our response, that may be relevant |54| . These 
effects arc currently under investigation and will be reported elsewhere. 

Our primary goal was to provide expressions for the differential cross sections that could be used in simulations 
of the early thermal evolution of neutron stars born in the aftermath of a supernova explosion. The microphysics 
of neutrino scattering affects the rate of diffusion, which in turn affects macroscopic observables such as the cooling 
rate and the neutrino emission from core-collapse supernovae. Our results, which extends both to the low and high 
temperature regions, are well suited for use in simulations of core-collapse supernovae and early thermal evolution of 
neutron stars. 
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APPENDIX A: CALCULATION OF THE POLARIZATION TENSOR FOR THE RELATIVISTIC 

SUPERFLUID 

In this appendix we detail the calculation of the medium polarization tensor for the one-component, relativistic 
superfluid. In Section A.l we evaluate the contribution from the quasi-free response function. In Section A. 2 we 
include the contribution from the collective response. 

1. Quasi- free response 

The contribution to the medium polarization tensor can be associated with the first diagram on the right-hand side 
of Fig. ^] That diagram has the value 

= -4 / ^Tr[r£S(P)r^(P + Q)]. (Al) 

The factor of \ is due to the doubling of fermion degrees of freedom in the Nambu-Gor'kov formalism. Evaluating 
the trace over Nambu-Gor'kov indices gives 

II£(Q) = -i\ J -^Tr[r(c v -c Al5 )G + (P)^(c v -c Al5 )G + (P + Q) 

-r(c v - c Al5 )E-(P)^(c v + c Al5 )E+(P + Q) 
-r(cv + c Al5 )E+(P)^(c v - c Al5 )E-(P + Q) 
+r(cv + c Al5 )G-{P)Y{cv + c Al5 )G-{P + Q)], (A2) 
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where the remaining trace is over Dirac indices. We will write down the value of each of the four traces above. First, 
Tr[y(c y - c Al5 )G+(Ph»(c v - c Al5 )G+(P + Q)} 

= b^-g][(p + g )g-g +q ] ^ + ^ )T + (P'P + q)- 2 ^>^ ( P ,P + q)}, 

where 7^ " is defined by 

(Pi, Pa) = Tr[ 7o yA± 7o r A± ], (A3) 

and Wjl" is defined by 

Wf (Pi,p 2 ) - ^[707^707^+75]. (A4) 

Explicit evaluation of T± v and W± K gives Eqs. 1)391 140(1 . respectively. The second term in l(A2(l is 

Tr[- 7 ^(cy - c Al r )E-(P)Y{cv + c Al5 )E+{P + Q)\ 
-A 2 



{(c 2 v + (? A )U1 V (p, p + q) + 2 CyCA Vf (p, p + q)} , 



^-^][(P+«)g-^+c,] 

where is defined by 

Wr(Pi,P2) = Tr[ 7 ^ i7 ^ 2 ], (A5) 

and V± is defined by 

Vr(Pi.Pa) - Tr[ 7 o7 /i A± 707^+75]. (A6) 
Explicit evaulation of Uj^ and V± gives Eqs. 1(4111421) . respectively. The third term in IA2I is 
Tr[-7^(cy + c Al5 )E+(P)Y(c v - c Alh )ET (P + Q)\ 

~ A {(4 + 4)^(p,p + q) - 2c V C A V fw (p,P + q)} • 



And the final term in (|A2|) is 

Tr[y(c y + c Al5 )G~(P)Y(cv + c Al5 )G-(P + Q)] 

(P0 — E p )(po + qo — E p+ q) f 2 2 . ™. . ..,»!/, . s-i 

= ^-g][(p + g )g-g +q ] ^ + C ^ T - CP,P + q) + 2«VCxWi: ( P ,P + q)} 



Next, we consider the Matsubara sums: 



A{b ,£,) - — y T — ■ — -2 — -P2TT2 — 7m y AJ ) 



and 



B(£ , ,.E) = 4S f7 I 1 , 2 u 2 772T ' ( A8 ) 

P „ [(Po + go) 2 - C J LPS - £' 2 ] 

where po — —i(2n + 1)tt/(3 is a fermionic Matsubara frequency; qo — — i2m7r//3 is a bosonic Matsubara frequency; and 



£' = v 7 ^' 2 + A 2 and £ = V# 2 + A 2 . 
Explicit expressions for A and -B are given in Eqs. 1(371 138|) . respectively. Adding all the terms yields Eq. ((36(1 . 
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2. Collective response 



In this section we detail the evaluation of the quantity depicted on the right-hand side of Fig. |H| corresponding to 
the collective response of the medium. First we evaluate the numerator of that quantity, then the denominator. 



a. The numerator 



The numerator of the right-hand side of Fig. [5] contributes 

~ G2 \\ J ff 4 ( 2^ 4 Tr[r£S(P)r g S(P + Q)} Tr[T H S(K)T z S(K + Q)\ = 

where 

J"(Q) = J ^Tr[T z S(P)T H S(P + Q)}. (A9) 

To evaluate P(Q) we first compute 

Ty[T z S(P)T h S(P + Q)} 




r(c v -c Al5 ) \ / G+(P) S"(P) \ 

-r(c v +c Al5 ) J \E+(P) G-(P) J 



vys\( G+(P + Q) ~-(P + Q) \ 
z 75 j \ E+(P + Q) G-(P + Q) J 

= Tiiir(c v -c Al5 )G + (P) l5 Z + (P + Q) 
+irf{cv - c Al5 )~- (P) l5 G+(P + Q) 
-iy(cv + c Al5 )E+(P) l5 G-(P + Q) 
-<y(cv + c j4 75)G-(P)7 5 S-(P + Q)]. 



Now, 



Next, 



And 



Finally, 



Tr[i 7 "(c y - C4 7 5)G+(P)7 5 S + (P + Q)] 
( Po +E p )(-iA) 



CvWf '(p, p + q) + caV^ u (p, P + q 



,,o. 



lt[iy*(<V - c A 7 5 )S-(P) 75 G+(P + Q)] 
(iA)((p + g) + P p+q 



b>g-g][(p + g)8-g +q ] 



cvW^°(p, P + q) + caV^ u (p, p + q 



,,o. 



Tr[-i 7 ^(cy + c^75)S + (P) 75 G-(P + Q)] 
(iA)((p + g)o-^ P+q 



l^-^][(P + 9)g-^ q ] 



CyZ^°(p, p + q) c A V^ U (p, p + q) 



Tr[-^(c y + c Al5 )G~ (P) 75 S- (P + Q)] 
( Po -E p )(-iA) 



b%-%][(p+Q)* -tM 



Cy^ (p,p + q)- CA V^ U (p,p + q 



no, 
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To proceed we evaluate the Matsubara sums 



r (F i F) = 1 y pq + e ' 

1 ' ' ~ P^Kpo+qoy-eM-? 2 } 



1 l\ (t\ (e'W ( Z-E 1 ? + E' \ 



1 »^-'«l^?e + 5^fe)> 



and 



C (F 1 F\ = 1 V" Po + qo + E 

[ ' j " /?^[(po + <Zo) 2 -»§-£' 2 ] 

y^7){[n(0-n« )] |^ 



+[1 _„ (s) _„ (f)1 (__!+£_ __!_£_), 



If we define 



C(E',E) = -C 1 (E',E) + C 2 (E\E) 



i /r ,~ e+c E-E'+z-e 

— {[n(0 -n(£ )] 



then an expression for 7 M is 



[c(E p ,E p+(l )[c V Uf(p, p + q) + c A Vf (p, p + q)] 
+ C(-E p , -£ p+q )M^°(p, p + q) - c A V?(p, p + q)]} . 
After performing the J" d 3 p integral, the terms proportional to ca vanish, leaving 

I"(Q) = J^^Ac v [c(E p ,E p+cl )Uf (p, p + q) + C(-£ p , -£ p+q ) Z^°(p,P + q) 
This is Eq. (gH. 
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b. The denominator 



In this section we evaluate the denominator of the right-hand side of Fig. H3 It has the value 
g(Q) = 1 + ^2^ ^Tr[S(P + Q)T H S(P)T H ] 

= 1 + G2 IJ 0^ns(P + Q)r H s(P)r H ] 



^,r(^Tr\ G+{ ■ p + Q ^ s-(P + Q)W o i75 \ 

J (2ir)* \Z+{P + Q) G-(P + Q)J O) 

( G+(P) S-(P)\ / i 7s \ 
\Z+(P) G-(P) )\i<y 5 J i 

= 1-gJ ^Tr[S-(P + Q) 75 S-(P) 75 + G + (P + Q)7 5 G-(P)75 
+G-(P + Q) 75 G + (P)75 + S+(P + Q) 75 S+(P) 75 ] 

= 1 + G / (2^^ + P-^ 

„ ((? + 9)0 + ^p + q)(P0 - Pp) - A" + ((p + g) - ^p + q)(j50 + E p ) - A 2 



1 1 



2\ 



> 90 + £p - £p+q <?0 ~ £p + £p+q 

[1 - n(Cp) - n(e P+q )] (e P +qe P + £p+q^p + a 2 ; 

1 1 



v QO - £p - £p+q 90 + £p + £p+q 

This is Eq. gBJ). 



APPENDIX B: COLOR-FLAVOR TRACES IN THE POLARIZATION TENSOR FOR CFL 



In Appendix A we evaluated the polarization tensor for the one-component, relativistic superfluid. This involved 
traces over Nambu-Gor'kov and Dirac indices, as well as a sum over Matsubara frequencies. Similar calculations 
are required to evaluate the polarization tensor for CFL. The two gap color-flavor structure of the quark propagator 
complicates things. We discuss in this appendix the color-flavor traces that must be computed to evaluate the CFL 
polarization tensor. First we discuss the color-flavor traces involved in the evaluation of the quasi-free response, and 
then we discuss the color-flavor traces involved in the evaluation of the collective response. 



1. Quasi- free response 



The contribution to the polarization tensor from the quasi-free response of the medium comes from the first term 
on the right-hand side of Fig. |H| Evaluating this diagram involves a trace over Nambu-Gor'kov and Dirac indices, as 
well as a Matsubara sum - these calculations are similar to those in the relativistic case. The complication in the 
present case comes from the trace over color-flavor indices, due to the two gap structure of CFL. The color-flavor 
trace involves the quantities 

P« = Tr[cyP m cyP n ] + Tr[c A P 
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] + Tr[ciP m cvP„], 
Rgl = Tr[cyP m CyP„] + Tr[c A P 

and 

R<$ n = Tr[c v P m c A P n ] + Tr[c A 
Explicit values for these quantities are given in Eq. (|55|l . 

2. Collective response 

We now want to calculate the contribution to the quark polarization tensor from the collective response associated 
with the H boson, depicted in Fig. El The color-flavor traces involved in calculating this contribution are much simpler 
than those involved in calculating the quasi-free response, since the -H-quark coupling does not mix the singlet and 
octet quasi-quarks. To evaluate the numerator of Fig. we only need the coefficients 

B$ = Tr[c y P m AfP m ] = Tr[c y P m MP m ] = Tr[c A P m MP m ] = Tr[c A P m MP la ], 

This is explicitly evaluated in Eq. Q59JI . To evaluate the denominator, we only need the coefficients 

R$ = Tr[MP m MP m ] = Tr[MP m MP m ]. 

This is explicitly evaluated in Eq. (|62|l . 
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